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Abstract 

Practical implementations of quantum computing are always done in the presence of decoherence. 
Geometric phase is useful in the context of quantum computing as a tool to achieve fault tolerance. 
Recent experimental progresses on coherent control of single electron have suggested that electron 
in quantum dot systems is promising candidate of qubit in future quantum information processing 
devices. In this paper, by considering a feasible quantum dot model, we calculate the geometric 
phase of the quantum dot system in nonuitary evolution and investigate the effect of environment 
parameters on the phase value. 
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The quantal geometric phase was first discovered by Berry |l| in 1984 in considering 
the quantum systems under cyclic adiabatic evolution. It has aroused much attention of 
researchers due to its importance. Since then the original notion of Berry phase has been 
extended to a general concept of geometric phase for pure states as well as for mixed states. 
The extension to pure states in nonadiabatic cyclic evolution was developed by Aharonov 
and Anandan Q in 1987, and that to pure states in nonadiabatic and noncyclic evolution 
was done by Samuel and Bhandari [s| in 1988. Further generalizations and refinements, 
by relaxing the constrains of adiabaticity, unitarity, and cyclicity of the evolution, have 
since been carried out 0, 0]. While all these extensions are of quantum systems in pure 
states, Uhlmann [6| was the first to address the geometric phases of mixed states within 
the mathematical context of purification. A physical definition of geometric phases for 
mixed states in unitary evolution was put forward by Sjoqvist et al. [7| in 2000 based on 



quantum interferometry, and it was recast in a kinematic description by Singh et al. 
The generalization of mixed geometric phases to quantum systems in nonunitary evolution 
was given by Tong et al. {9} in 2004. More works on geometric phases related to states for 
open systems may be seen in Refs. 101 ] . 

The geometric property of the geometric phase has stimulated many applications. It has 
been found that the geometric phase plays important roles in quantum phase transition, 
quantum information processing, etc. The geometric phase shift can be fault tolerant 
with respect to certain types of errors, thus several proposals using NMR, laser trapped 
ions, etc. have been given to use geometric phase to construct fault-tolerant quantum 
information processer 12], and the fault-tolerant geometric quantum computation gate has 
been demonstrated in experiments using NMR 13]. 

Geometric phase is useful in the context of quantum computing as a tool to achieve 
fault tolerance. Practical implementations of quantum computing are always done in the 
presence of decoherence. Fortunately, recent experimental progresses on coherent control of 
single electron have suggested that electron in quantum dot systems is a promising candidate 
of qubit in future quantum information processing devices [14J, because it has long spin 
coherence time. This start us to investigate the geometric phase of quantum dot systems in 
nonuitary evolution. In this paper, we calculate the geometric phase of a feasible quantum 
dot model and investigate the effects of the environment parameters to the phase value. 

The model is illustrated as FigJTJ Two quantum dots, QDi and QD 2 , are coupled to each 
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FIG. 1: Illustration of the model. 

other with strength s\. An electron is trapped in the quantum dots and it tunnels between 
the two quantum dots. Only one energy level is considered in each quantum dot, and hence 
the electron and the two quantum dots construct a two-level quantum system, a qubit. The 
environment of the system consists of another quantum dot, QDo, and two leads connecting 
to QD . The left lead has higher chemical potential than the right lead. Electrons can 
tunnel from the left lead to QD and then tunnel out to the right lead. For simplicity, we 
assume that only one electronic state with the energy level Eq in QDo is correlated and 
fiL > E > fi R , where and are the chemical potentials of the left lead and the right 
lead respectively. Once there is an electron in QD , it will affect the coupling between QDi 
and QD2 by changing the coupling strengthes from s\ to S2- This is an interesting model, 
of which the relaxation and decoherence and quantum measurement have been well studied 
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161 ]. The model is easily performed in experiment, and it may play a potential selection 
for geometric quantum computation of using quantum dot systems. 

Noting that the qubit system, comprising the trapped electron and the two dots, is an 
open system being in mixed state, we use the formula of geometric phases for mixed states 
in nonunitary evolution given in Ref. Q. For an open quantum system, described by the 
reduced density operator, p(t) = Y^t=i u k{t)\(j)k{t)){4>k{t)\, t G [0, r] the geometric phase is 
given by the formula, 

7(r) = Arg[^v/^(O)^(r)(0 fc (O)|0 fc (r))e-/o T ^WI^(^A (1) 
^ k=i ' 

where Ukit) is the k— th eigenvalue of the reduced density matrix, \4>k{t)) is the corresponding 
eigenvector, and r is the total evolutional time. 

In order to calculate the geometric phase of the qubit system, we need to obtain the 
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reduced density operator. The Hamiltonian of the large system can be expressed as 



H =H S + H e + Hi, 

H s =E\a\a\ + E 2 a\a 2 + S\{a\a 2 + a 2 ai), 

H e =E clc + 2_. Eic\ci + 2j E r S r C r + /^(^V]^ + + ^ C ) 



Lr 



Hi = (s 2 - si)cj,c (ala 2 + a\ai). (2) 

Here, H s , H e , Hi are the Hamiltonians corresponding to the system itself, the environment 
and the interaction between the system and its environment, respectively; a\ and a\ (a± and 
a 2 ) are the electron creation (annihilation) operators in the two quantum dots; c\ and c\. (q 
and c r ) are the electron creation (annihilation) operators in the environment corresponding 
to the left lead and the right lead respectively; E\ and E 2 are the energy level of QDi and 
QD 2 ; (fi r ) is the coupling parameter of left (right) lead with the quantum dot QD . For 
simplicity, we have considered electrons as spinless fermions, and we have used Ei, E r , Qi, 
and Q r to represent E Lh E Rr , VL Lh and Q Rr respectively. 



The wave function of the large system, |^/(t)), satisfies the Schrodinger equation, i 



dm)) 

dt 



H(t)\^/(t)) . The reduced density operator p(t) may be expressed as the partial traces of 
| l I / (t))(\[ / (t)| with respect to the environment consisting of the quantum dot QD and the 
two leads, p{t) = tr Do g(t), where g(t) = tr^ s | 1 i r (t))( 1 i r (t)|. Following the method used in 



Refs. 16|, we may get the equations of motion for the elements of density matrix g(t). The 
bases of g(t) consists of four discrete states, |1) = |1,0,0), |2) = |1,0, 1), |3) = |0, 1,0), 
|4) = |0, 1,1), where \rii, n 2 , n 3 ) means that there are rij, n 2 , n 3 electrons in QD 1; QD 2 , 
QD respectively. In the approximation of constant density of states, let Tl = 2tc\Ql\ 2 pl 
and T^j = 2^1^^^^, where pl (pr) is the density of states for the left (right) lead, and Ql 
(fin) denotes the constant coupling parameter fl[ (fl r ). T L (T R ) depicts the tunneling rate 
between the left (right) lead and QD . In this case, the elements Qij of the density matrix 
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g(t) satisfy [171] . 
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The initial condition is taken as Qij\t=o = 1, for % = j = 1, or 0, for all other i, j, correspond- 
ing to the case that the electron is in QDi and no electron is in QDo- Here eo = Ex — E 2 , is 
the energy difference of the energy levels of QDi and QD2. The elements pij of the reduced 
density matrix p(t) of qubit can be then expressed as 



Pll = 1 - P22 = 011 + 022, Pl2 = P2I = 013 + 024- 



(4) 



Once the reduced density matrix is obtained, we can calculate its eigenvalues U)k(t) and 
eigenvectors \<f>k(t)), and we have 

1 ± a/ (pxx - P22) 2 + 4|pi 2 | 2 
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The initial condition taken above implies ux(0) = 1, uj 2 (0) = 0, and |0i(O)) 
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The evolution of the system can be illustrated by the path traced in Bloch sphere. The 
three-dimensional coordinates in the Bloch sphere are x = p 12 + P21, V = i(pn — P21), 
and z — pxx — P22, respectively. By using the Four-order Runge-Kutta method, we may 
numerically resolve the differential equations in ([3]) and obtain the value of the density 
operator. Fig. [2] shows the path traced by the state of the system, where the parameters 
are chosen as Tl = 1.0, Tr = 2.0, sx = 1.0, s 2 = 0.5, eo = —2.0. Hereafter, we take 



the parameter si as the base unit. All the other parameters with energy dimension, such as 
Tl, r#, s 2 , are measured by the unit S\, and the time is measured by l/s\. As the time goes 
on, the path starts from (0,0, 1), which corresponds to the state that the trapped electron 
is in QD 1; and moves spirally to (0,0,0), which corresponds to the state that the electron 
has half probability in QDx and half in QD 2 . 




FIG. 2: The Bloch sphere of the density matrix. 

Substituting Eq. (jSJ) into Eq. ([ID, we can calculate the geometric phase of the system. It 
may be simply expressed as j(t) = i ((pi(t)\4>i(t))dt. To sketch out the changing tendency 
of the geometric phase, we numerically calculate the geometric phase. The parameters are 
again chosen as T L = 1.0, T R = 2.0, Si = 1.0, s 2 = 0.5, e = —2.0. The result is shown 
as Fig. El The geometric phase is usually put in region [0, 2tt) (mod 2n). In order to show 
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time (units of 1/sJ 

FIG. 3: The geometric phase as a function of time. 

entirely the changing tendency of the phase and express clearly the path dependence of the 
geometric phase, here we give the schematic by using the calculated values without making 
a 27r-modulus. The recast of the results in [0, 2tt) is trivial. 

From Fig. [31 we find that the geometric phase is changing as the time is going on, and 
it finally saturates to a constant value. The saturation value is a characteristic value for a 
given configuration of parameters, which may be simply called as the characteristic geometric 
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phase (CGP). This is consistent with the 'geometricity' of the geometric phase, that is, the 
geometric phase is only dependent on the path traced by the state of the system, but not on 
the dynamics. When evolutional time is small, the spiral path has large spiral radius and 
the changing of the path is notable, and thus the changing of the geometric phase is obvious 
too. With evolutional time going on, the spiral radius of the path becomes small and the 
changing rate of the spiral path are reduced, and therefore the changing of the geometric 
phase will be reduced too. The system will finally evolves to the point (0,0,0), and from 
then on the path will be little changing, and so does the geometric phase. 

In the model, there are three environment parameters S2, Tl, Fr- We now investigate 
the effects of these parameters on the phase values. For this, we will consider two kinds of 
geometric phase values, the geometric phase corresponding to the whole evolutional time, 
i.e., the CGP, and the geometric phase corresponding to a special time interval T. 

Firstly, we observe the effect of the parameters on the CGP. Fig. [U^a) shows the effect 




FIG. 4: The geometric phases, CGP and j(T), as functions of the parameters S2, Tl, and Tr. The 
parameters except for the one taken as variable are chosen as Tl = 1.0, Tr = 2.0, s% = 1.0, S2 = 
0.5, e = -2.0. 

of S2 on CGP. From the figure, we see that CGP is strongly dependent on the parameter 
S2- Specially, CGP is infinitely large at S2 — si. This is a reasonable result, because S2 = s± 
means that the environment does not affect the qubit system. In the case, the qubit is in 
the pure state, which is evolving repeatedly along a closed circle in the Bloch sphere, and 
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CGP will accumulate infinitely as the time is going on. However, as the parameter s 2 — s% is 
becoming large from zero, the value of the phase will reduce. The phase values will approach 
to zero when s 2 — Si is large enough. This may be explained by the following argument. The 
larger s 2 — Si means the larger correlation between the environment and the qubit system, 
which leads to the smaller spiral radius of the path traced by the state of the system. When 
the environments' effect is stronger enough, the path may approaches to a line directly from 
(0, 0, 1) to (0, 0, 0) and the corresponding geometric phase will be near to zero. 

Figs. H] (b) and (c) show the effect of parameters I\ and Fr on CGP. From the figures, we 
find that the two curves in the figures are similar. CGP becomes infinitely large at Tl = 
or Tr = 0, and it is also approaching to infinity as Fl or Fr is going to large values. These 
observations are consistent with the physical construction in the model, as we have taken 
Fr : = 2 in Fig. 0] (b) and Fl = 1 in H] (c). Roughly speaking, when is small and Fr 
is large, electrons are hard to tunnel into QD from the left lead but easy to tunnel out of 
QDo- There is nearly no electron staying in QDo in all the time, i.e., the coupling between 
QDi and QD2 is mainly s±. The effect of the environment on the qubit is negligible, and 
the qubit may be taken as a closed two-level system with coupling strength si. The picture 
of CGP corresponding to the case is the left part of Fig. HKb) or the right part of 0(c). 
When is large and Fr is small, electrons are easy to tunnel into QDo from the left lead 
but hard to tunnel out of QD . There is an electron staying in QD almost in all the time, 
i.e., the coupling between QDi and QD2 is dominated by s 2 . The effect of the environment 
on the qubit is only to change the coupling strength between QDi and QD 2 from si to s 2 , 
and the qubit system may be taken as a closed system but with coupling s 2 . The picture 
corresponding to this case is the right part of the curve in Fig. H](b) or left part of the curve 
in Fig. H](c). When Fl and Fr are in the same order, the qubit is an open system in mixed 
state. The path traced by the mixed state is a spiral curve and so corresponds to finite 
values of CGP. 

Secondly, we observe the effect of the parameters on the geometric phase for the special 
time interval T. If there is no coupling between the qubit and the environment, or s 2 = s\, 
the qubit system will be in a pure state and it will evolve from the initial state (0, 0, 1) back 
to itself after a time interval T, making up a closed circle in the Bloch sphere. In the case 
where eo = —2 and s\ = 1, we have T = n/\/2, and the geometric phase corresponding to 
the closed cycle is 7(T) = n — n/\/~2. However, if s 2 7^ si, the path traced by the state in the 
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Bloch sphere will become an unclosed curve and the geometric phase j(T) will be changed 
under the effect of the environment. Therefore, 7(T) may be used to describe the effect of 
the environment on geometric phase in an finite time, during which the pure state evolves 
one circle. Figs. 11(d), H|e) and|U(f) show the effect of parameters S2, and Tr on 7(T), 
respectively. The curves in the figures may be explained by applying a similar discussion as 
above. 

In conclusion, we have calculated the geometric phase of a feasible quantum dot model and 
investigate the effects of the environment parameters to the phase value. Here, we not only 
presented the parameters's effect on the characteristic geometric phase, which corresponding 
to the whole evolutional time, but also studied their effect on the geometric phase in a finite 
time interval T, defined by using pure state without the effect of environment. The approach 
of calculating the geometric phase in the paper is reliable. While the other approaches 
of defining the geometric phase of open systems have met criticisms [lfl, the kinematic 
approach used in the pap er has been widely applied to investigate the open systems in 
various environments 19J. Our investigation on geometric phase is helpful to completely 
understand the properties of the quantum dot system. 

This work is supported by NSF of China under Grant Nos. 10675076, 10875072 and 
10804062. 
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